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Abstract 

We explore some of the gravitational features of a uniform in- 
finitesimal ring both in the Newtonian potential theory and in Gen- 
eral Relativity. We use a spacetime associated to a Weyl static so- 
lution of the vacuum Einstein's equations with ring like singularity. 
The Newtonian motion for a test particle in the gravitational field 
of the ring is studied and compared with the corresponding geodesic 
motion in the given spacetime. We have found a relativistic pecuHar 
attraction: free faUing particle geodesies are lead to the inner rim 
but never hit the ring. 

PACS numbers: 02.40.-k, 04.20.-q, 04.20.Jb, 04.40.-b, 04.25.-g 

1 Introduction 

The astrophysical and elementary particle physics importance of ring Hke con- 
figurations is evident: there are several ring like structures as for example in 
galaxies P] and planets; and the lowest energy state of a closed string is a ring. 
However, there are only few solutions of the Einstein's field equation that rep- 
resent the gravitational field of a ring. With somehow constrained hypothesis 
of static configuration we mention the Weyl-Bach solution which is the general 
relativistic analog of a Newtonian ring of constant density [2| , given in terms of 
elliptic functions. The ring is not a simple line source P] as it will be clear from 
the strange effect they have on the particle motion. The solution is asymptoti- 
cally fiat but the outer communication region is not simply connected. 

We should mention several studies of self-gravitating Newtonian rotating 
rings. The rotation is of primordial importance to the ring's stability (and 
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instabilities as well) but it has no effect on the Newtonian potential and even in 
General Relativity its (magnetic part of the curvature) contribution to geodesies 
is usually weaker than the sole static (electric part of the curvature) one. 

The purpose of this paper is to study some properties of the static ring solu- 
tion due to Weyl-Bach j2j using geodesies. In particular we shall be interested 
in the attractive or apparent repulsive character of the ring singularity as well as 
its directional feature. Static axially symmetric solutions of the Einstein's equa- 
tions usuahy are characterized by the presence of string hke singularities (conic 
singularities) and its higher dimensions generahzations 0. In some cases, these 
singularities arise as supporting devices of an otherwise dynamical configuration 
of masses. 

The ring has a different attraction at its rim. Free falling particles are led 
to the inside part of the Weyl-Bach ring. This might be linked to the different 
tension needed to keep the ring static. Besides that the Weyl coordinates usually 
compact singularities and event horizon into a lower dimensional region so that 
the physical distances are not so evident. Of course the ring is black since the 
lapse function vanishes at the ring and the spacetime is static. But the curvature 
invariants diverge there as well. 

In Sect. |2lwe present the equations to be solved. In subsect. 12. II we present 
a summary of the main expressions associated with the static axially symmetric 
spacetimes solutions of the Einstein's field equations and the geodesic equations 
for test particles evolving in these spacetimes. Then, subsect. 12.21 shows the 
associated Newtonian equation of motion for comparison. At last, expected 
behavior for axial and plane motions is given in subsect. 12.31 

Next in Sect. Elwe show the potential associated to the ring. We paid special 
attention to the motion of test particles: The Newtonian motion (section^ll and 
geodesic motion (sectional in the spacetimes associated to the Weyl solution 
are considered. Most of the motion is not trivial. 

We think it is important to understand the effect of the ring itself on the 
particle motion to appreciate the more complex configurations in which a ring 
is just part. 

2 The Dynamical System of Equations 

Let us fix our coordinate system with axial symmetry: let r > be the coordi- 
nate away from the axis and z € M the coordinate along the axis. The dynamical 
system evolves with either a time coordinate or a proper time. In any case a 
curve in the half-plane r—z will be parametrized by a "time" parameter and its 
coordinate rate of change denoted by r, i and so on. 

The dynamical system of equations wih come either from the geodesic equa- 
tions in a spacetime or from Newton equations of motion of a test particle in 
gravitational potential. We assume that both the spacetime and the potential 
are static and axisymmetric. So our dynamical system depends either on metric 
functions for geodesies or on the gravitational potential for Newtonian motion. 

Next we present the exphcit equations to be solved. 
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2.1 Geodesies in Weyl Solutions 



The static spacetime of an axially symmetric body can be described by the Weyl 
metric 



-20 r^27 



{dr 



d7j 



■dip' 



(1) 



where the functions (p and 7 depend only on r and z; the ranges of the coordi- 
nates (r, z, ip) are the usual for cylindrical coordinates and f £ R. The vacuum 
Einstein's equations {Rfi„ = 0) reduce to the Laplace equation in cylindrical 
coordinates, 

(p^rr + -4>,r + 4',zz = (2) 



and the quadrature, 



dj[4>] = r [((/. 2 _ dr + 20,,0, dz] 



(3) 



If 4> satisfies the Laplace equation Q then 7 is twice differentiable. The function 
4> determines the Weyl solution uniquely up to a constant. 

The geodesic equations in Weyl spacetimes have two constants of motion 
associated to the cyclic variables t and ip, 



E 



(4) 



Now the overdots mean derivation with respect to the proper time, t — s. The 
other two second order evolution equations are 



z = (f^ - i^)(7,2 - (t),z) - 2ri(7^r - 4>,r) - 4',zQ, 



(5) 
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Prom the constant of motion gfj.i, x^i" = 1 we find 
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This function has some similarity with the Newtonian effective potential. The 
motion is allowed only where F{r, z) < 0. 
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2.2 Newtonian Motion in Axisymmetric Potential 

In Newtonian gravitation, the motion of a test particle in a field of forces de- 
scribed by an axially symmetric potential (/), solution of Laplace equation, is 
characterized by two constants of motion: the energy H and the angular mo- 
mentum L, 

L = r^^, (8) 

where the over dots mean derivation with respect to the Newtonian time; and 
two second order differential equations 

f= — z = (9) 
From these equations we have the constant of motion 

H=]^(i^ +P) + V{r,z), 
where V{r, z) is the effective potential, 

V{r,z)^^ + ct^ir.z). (10) 
In the Newtonian case the motion is allowed where V{r,z) < H. 



2.3 Axial and Equatorial Motions 

It is easy to see that if the problem has both axial and planar symmetry at 
z = and no source except for the ring, then a test particle with z = = i is 
confined in the plane z = since the partial derivatives of the functions along 
z at z = is zero. The so called equatorial motion. 

There is also an axial motion. A test particle with r = = r has vanishing 
angular momentum L = and by hypothesis of axial symmetry and no source 
at the axis, the functions has vanishing partial derivative along r at r = 0. 
Therefore r — and the particle stays on the axis. 

Of course the origin is an equilibrium point as long as the hypothesis above 
are satisfied. 

The general motion of a test particle can be very compHcated for non trivial 
functions. Care must be taken because the functions have singularities at the 
ring. 



3 Ring Potential 

The ring we are concerned with uses the function (fi which solves the Laplacian 
everywhere in the half plane r — z except for the ring at z = 0, r = a. It is 
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Figure 1: Contour plot of the potential (j) for a constant mass density ring with 
mass M = 1 and radius a = 1. 
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the gravitational potential itself for a Newtonian motion or the metric function 
which sets the Weyl solution of the spacetime in which we have to solve the 
geodesic motion. 

The Weyl-Bach solution has as Newtonian image the usual potential for a 
ring of uniform density which is a solution of Laplace equation. It can be written 
as 

, K fMl) . f- .... ../^ . , (U) 



irRa \ Ra J 27r 7o ^7-2 + z2 + a2 - 2ar cosi^ 

where — {a + rY + and K(a;) is the complete elliptic integral of the first 
kind for x G [0, 1). The ring is located on the plane z = and its center has 
coordinates r — z — ^. We shall take from now on M = 1 = a. The contour 
plot of this function is depicted in Fig. 

The evaluation of the elliptic integral as well as its derivatives was made 
with an algorithm adapted from the one presented in 22l- Of special interest 
are the values of at the axis r = and at the plane z = 0. Using the 
integral representation above we get 



r^O 27r ./n \y2 I „2l3/2 



and 



lim (6 y = — / lim 




*o ' 27r 7o [r2 + + - 2ar cos (^]^/^ 

Thus, the only Newtonian source for the potential is the ring itself 
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Figure 2: Effective potential of the ring at z 
z = 0, i, 1. 



= 0, r = 1, for r G (0,3) and 




4 Newtonian Motion 

The equiHbrium (circular) position of a test particle, if it exists, obeys the 
equations f = i = f = z = 0, that is, from the effective potential ifTflll 

= , dz<i> = 0. 

The planar symmetry implies the position may be in the z = plane. And for 
the ring potential c^z^ > for z 7^ 0. Thus the motion is stable about the ring's 
plane for appropriate energies. 

In the disk inside the ring we have < so an equilibrium point is possible 
only for L = at r = but it is not stable since any amount of angular moment 
L will push the particle towards the inner part of the ring. 

If there is a velocity in the z directions, i 0, and for small values of angular 
momentum, the test particle moves up and down the disk inside the ring and 
for slightly higher angular momenta it may cross in and out of the ring. 

For L ^ 0, the centrifugal force pushes the particles away from the axis. 

Outside the ring but in its plane, there is a lower bound of angular mo- 
mentum, let us say Ls, beyond which there are stable equilibrium motions at 
a distance > a. Using and a = 1 = M we obtain Ts and Lg by setting 
dVldr = = d'^Yldr'^ at z = 0. We find that = 3.8396 and r, = 1.6095. 

We show in Figs. |2ti,b,c the effective potential for the ring M = 1 = a, for 
different values of angular momentum, L = 2~''Ls, for /c = 3, 2, 1 respectively. 

The motion of a test particle initially at rest at r = z = 1 { L — 0) is study 
in Fig. El 

In Figs. 0^,b we present trajectories with initial conditions r — z — 1, 
f — z = and different angular momenta, L — 1 and L — 4.5, respectively. The 
orbits in these two cases are bounded in a tri-dimensional region of the space. 
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Figure 3: Trajectory of a free falling particle in the gravitational field of a ring 
of constant density. The particle initially at rest {L = 0) at {r,z) = (1, 1). 
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Figure 4: Trajectory of a free falhng particle with initial conditions r = z = 1, 
r = z = 0. 
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(a) Case L = 1. (b) Case L = 4.5. 
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Figure 5: Contour plot of 7 for the Weyl-Bach solution associated to a ring with 




5 Geodesies in Weyl-Bach Ring 

Now for the general relativistic case we have two "potentials", ■0 and 7. The 
metric functions for the Weyl-Bach ring are better expressed in toroidal coor- 
dinates (77, ^) € [0, 00) X [0, 2tt) which are related to the axial coordinates (r, z) 

by 

sinh rj sin ^ 



cosh ?7 — cos f coshT^ — cos^' 

^ z'^ — 0?" 
cot ^ — , coth 77 



7-2 + + 



2az 

The Newtonian potential can be cast as 



2ar 



(j) = —ae ''/^K(k) cosh 77 — cos ^, 

where K(k) is the complete elliptic integral of first kind, = 1 — e~'^^ and 
a = Then the 7 function for the Weyl-Bach solution is P ITTIITII . 



-K 



K M + k'- 



(2 



cos^ 
sinh 77 



- 2E U - K 



, cos^ 
sinh 77 



(12) 



where E is the complete elliptic integral of second kind and k! = e~^. For the 
actual computations of the function 7 this formula is not particularly interesting. 
We found more convenient its evaluation by direct integration of equation ^ 
using Gauss quadrature methods for chosen integration paths. The contour of 
the function Ijl2|l is shown in Fig. 

In Fig. we display time- like geodesies of test particles moving in the grav- 
itational field of Weyl-Bach ring. The initial conditions are (r, z) — (2, 0) with 
initial velocity only in the z direction determined by the values of E shown in 
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Figure 6: Geodesies in Weyl-Bach ring solution with initial conditions (r, z) = 
(2,0), r — and some values of E with a = I = a. 
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the graphic and L = 0. For low values of E, we have a clear repulsion at the 
beginning of the trajectories. The repulsion happens a little later for higher 
values of E. The geodesic with E = 3 seems to suffer only attraction. Therefore 
we have a separatrix closer to this value of E. 

These geodesies have interesting asymptotic behavior. They approach the 
ring from inside its interior at very high values of the proper time, that is, the 
Weyl-Bach ring has a directional singularity and the particle hits the ring only 
at infinite proper time. This behavior seems to be generic. 

The behavior of these geodesies is characterized by the following two features: 

a) All these geodesies are asymptotically 'radial'. That is, there is a privileged 
directions around the ring, pointing to a directional singularity of the Riemann 
tensor similar to the singularity of the Chazy-Curson metric jH]. The curvature 
scalar invariants wi = ^CabcdC'""^ and W2 = i^C'a&C'cd C'e"/' ^^ere C'"''^ is 
the Weyl tensor have different limits when one approaches the singularity from 
different directions 0. Both invariants blow without bound when the limit is 
taken from the interior of the ring and approach zero from up or down directions. 

b) We have the 'freezing' of the motion in Weyl coordinates. The numerical 
computation of these geodesies for very high values of proper time indicates 
that they will reach the singularity at an infinite proper time. 

In Fig. 13 we plot geodesies initially at {r, z) = (2,0), with E = 1, initial 
velocity with no component in the r direction, and L = 1,25,50,75, and 100. 
All the geodesies finish in the ring by its inner disk at infinite proper time. 
Finally, in Fig.|Hlwe present five trajectories of test particles initially at rest at 
(r, z) = (1, 1), (2, 1), (2, 0.01), and (2, -0.5). The geodesic that starts at (2, 0.01) 
suffers a very strong repulsion. Again, all the geodesies fall in the ring exactly 
as the previous case. 

It is clear the structure of the inner part of the ring is non trivial. Any 
geodesic approaching the inner part of the ring has either infinite proper time 
or infinite proper length. One can compute a proper distance to the ring of a 
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Figure 7: Geodesies in the Weyl-Bach ring solution starting at {r, z) = (2,0) 
with E = 1 and several values of L. 




Figure 8: Geodesies of Weyl-Bach ring solution starting at rest from (r, z) = 
(1,1), (2,1), (2.00,0.01), (2.0,-0.5). All the geodesies reaeh the ring singularity 
from the inner disk. 
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curve in the (r, z) plane, in which the ring is at (a, 0): 



(a,0) 



(13) 



Let us show that it diverges according to the different directions of the ap- 
proach using the toroidal coordinates toroidal coordinates. As ^ cjo 
and ^ € ^) the approach is from inside. Using the asymptotic of the EUip- 
tic functions one can show that 



lim (7 — 0) = — hm 

(r,z)— ►(a,0) 



— coshr/ cos 4 

flTT / 



Thus the ring is at finite distance when approached from outside (in which 
cos^ > 0) because the integrand goes to zero, whereas it is at infinite distance 
when approached from inside (in which cos^ < 0) because the integrand di- 
verges. Therefore a fixed physical distance to the ring means greater (smaller) 
coordinate distance in the directions where lim5]^o(7— 0) is small (large). Hence, 
the particles should appear "repelled" ("attracted") in the directions from where 
the ring is physically nearby (far away). This agrees with what can really be 
seen in the figures. 



6 Discussion 

We investigated the gravitation induced by a ring both in the Newtonian and 
in General Relativity dynamics of test particles. Although related, the space- 
time associated to the Newtonian potential of a ring has quite distinct features. 
We have learned that line sources in General Relativity exhibits directional 
singularities and the results above for the Weyl-Bach ring give explicit example 
of them. 

Furthermore, the Weyl coordinates have the tendency to compact a whole 
region into a singularity. This is the case for the Schwarzschild solution in which 
the event horizon and the physical singularity (with the topology of S*^ x R and 
R respectively) are displayed in Weyl conformal coordinates as a single world 
finite line segment (/ x M where / C M is a line segment). 

And we find very interesting that the inner side of the ring is very attractive 
but is not accessible for particle geodesies because it is too far away. The test 
particles approach the singularity in a privileged way: They arrive along radial 
directions of the ring inner disk. We presented geodesies which take an infinite 
amount of proper time to hit the ring. This happens also in the extreme case 
of the Reissner-Nordstrom solution. The physical distance to the event horizon 
is infinite. 

The geodesies display gravitational field with apparent repulsive regions. 
This can be either a coordinate effect as pointed out above or may indicate the 
presence of very high tensions in the ring. Probably both. 
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The apparition of tensions in the Weyl solutions is a known fact j^]. The 
imposition of a static geometry and the Einstein's equation creates some devices 
Hke strings, struts, membranes, etc. to support an otherwise dynamical config- 
uration. In the present case we have some kind of strong hoop tension along 
the ring. As we know the spacetime is sensitive to both density and pressure or 
tension while in the Newtonian gravitation the density of the source suffices for 
the gravitational potential. 

So far, there are few self-gravitating ring solutions of Einstein's equations 
[13] . The reader should be cautious about rings solutions in the literature (HI 
Some have misprints, others have misinterpretation (see Appendix A). 
Nevertheless they are very interesting. 

In this paper we show some interesting behavior of test particles about a 
ring alone. We think the understanding of the gravitation of the ring itself is 
useful for configurations in which the ring is an important part. 

We thank FAPESP and CNPq for financial support and Eduardo Gueron 
for the critical reading. 

A. Rings Problems in the Literature 

In 1887 a very simple solution of Poisson equation which is singular at a ring 
was pubhshed And indeed it was interpreted as a ring-like configuration of 
matter. This interpretation is present, at least, in two important books |15lll(i) . 
Based on that, Letelier and Oliveira \J[ developed a series of new potentials with 
singularities at the ring and were misinterpreted as generated by rings. Moreover 
they were able to get a Weyl class of new solutions for Einstein equations in 
vacuum with axial symmetry. Nevertheless, as Gleiser and Pulhn pHl correctly 
showed, Appell solution is not just a ring. Actually there is a surface mass 
density in the plane of the ring. Let us compute it because this is the source of 
another mistake caused by a misprint. 
The Appell potential is 



where r and z are standard cylindrical coordinates and the constants M and A 
may take complex values. Physical potential is the real part of (/>. We promptly 
see the singularity at r = a, z = if A = ia, that is, it is singular at the ring. 
Nevertheless the z derivative is not null at z = 0. Actually there is a jump 
across this plane: 



If one takes the real part of the potential with A — ia, where a and M are 
positive real constants one gets the surface mass density: 



M 



+ {z- A)2 



lim [0,^] = -2MA (r^ + A^) ^ . 




for r < a. 
for r > a. 
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In the power is misprinted as — ^. It happens that there is a known disk 
with such a surface mass density, in the class of the so called Morgan & Morgan 
family ^21 of disks in General Relativity. 

This mistake lead us into another misinterpreted result. If both disks had 
the same surface mass density one could subtract one from the other leaving 
just the ring [^. But Appell disk (sic) and Morgan & Morgan disk do not have 
the same surface mass density! 

Of course, looking back, one could not have more than one ring as solution 
of Laplace equation with axial symmetry. There exist several theorems prov- 
ing the existence and uniqueness of solution of the Laplace equation. For an 
infinitesimal ring with axial symmetry, there is no possibility of other solution 
but the constant linear mass density. And the Weyl solution linked to the po- 
tential is also unique. This is the ring of the main part of this paper with the 
corresponding Weyl Bach solution spacetime. 
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